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This  experimental  study  investigated  the  preignition  reactivity  behavior  of  a  gasoline  surrogate  in  a 
pressurized  flow  reactor  over  the  low  and  intermediate  temperature  regime  (600-800  K)  at  elevated 
pressure  (8  atm).  The  surrogate  mixture,  a  volumetric  blend  of  4.6%  1-pentene,  31.8%  toluene,  14.0% 
n-heptane,  and  49.6%  2,2,4-trimethyl-pentane  (iso-octane),  was  shown  to  reproduce  the  low  and 
intermediate  temperature  reactivity  of  full  boiling  range  fuels  in  a  previous  study.  Each  of  the  surrogate 
components  were  examined  individually  to  identify  the  major  intermediate  species  in  order  to  improve 
existing  kinetic  models,  where  appropriate,  and  to  provide  a  basis  for  examining  constituent  interactions 
in  the  surrogate  mixture.  n-Heptane  and  1-pentene  started  reacting  at  630  K  and  640  K,  respectively, 
and  both  fuels  exhibited  a  strong  negative  temperature  coefficient  (NTC)  behavior  starting  at  700  and 
710  K,  respectively.  Iso-octane  showed  a  small  level  of  reactivity  at  630  K  and  a  weak  NTC  behavior 
starting  at  665  K.  Neat  toluene  was  unreactive  at  these  temperatures.  The  surrogate  started  reacting  at 
630  K  and  exhibited  a  strong  NTC  behavior  starting  at  693  K.  The  extent  of  fuel  consumption  varied  for 
each  of  the  surrogate  constituents  and  was  related  to  their  general  autoignition  behavior.  Most  of  the 
intermediates  identified  during  the  surrogate  oxidation  were  species  observed  during  the  oxidation  of 
the  neat  constituents;  however,  the  surrogate  mixture  did  exhibit  a  significant  increase  in  intermediates 
associated  with  iso-octane  oxidation,  but  not  from  n-heptane.  While  neat  toluene  was  unreactive  at 
these  temperatures,  in  the  mixture  it  reacted  with  the  radical  pool  generated  by  the  other  surrogate 
components,  forming  benzaldehyde,  benzene,  phenol,  and  ethyl-benzene. 

The  observed  n-heptane,  iso-octane,  and  surrogate  oxidation  behavior  was  compared  to  predictions  using 
existing  kinetic  models.  The  n-heptane  model  reasonably  predicted  the  disappearance  of  the  fuel,  but 
overpredicted  the  formation  of  several  of  the  smaller  intermediates.  The  iso-octane  model  significantly 
overpredicted  the  reaction  of  the  fuel  and  formation  of  the  intermediates.  The  1-pentene  model 
reasonably  predicted  the  fuel  consumption,  but  underestimated  the  importance  of  radical  addition  to 
the  double  bond.  The  results  of  this  study  provide  a  critical  experimental  foundation  for  the  investigation 
of  surrogate  mixtures  and  for  validation  of  kinetic  models. 

©  2008  The  Combustion  Institute.  Published  by  Elsevier  Inc.  All  rights  reserved. 


1.  Introduction 

Understanding  the  chemical  processes  that  cause  a  fuel  to  au- 
toignite  is  critical  for  solving  the  problem  of  knock  in  spark  igni¬ 
tion  engines  and  for  the  development  of  Low  Temperature  Com¬ 
bustion  engine  technologies  such  as  Homogeneous  Charge  Com¬ 
pression  Ignition  (HCCI)  and  its  derivatives.  As  the  autoignition 
phenomena  is  strongly  dependent  on  the  oxidation  chemistry  of 
the  hydrocarbon,  which  depends  on  the  hydrocarbon  structure 
and  the  temperature  and  pressure  of  the  combustion  environment, 


*  Corresponding  author.  Fax:  +1  (215)  895  1478. 

E-mail  address:  dmiller@coe.drexel.edu  (D.L.  Miller). 

1  Current  address:  Praxair  Inc.,  175  East  Park  Dr,  Tonawanda,  NY  14150,  USA. 


a  fundamental  understanding  of  the  low  and  intermediate  tem¬ 
perature  chemistry  of  the  different  hydrocarbon  structural  groups 
present  in  automotive  fuels  is  necessary.  With  the  increased  in¬ 
terest  in  the  effect  of  fuel  formulations  on  controlling  engine  op¬ 
erations,  extensive  studies  of  the  combustion  of  fuel-size  hydro¬ 
carbons,  neat  and  in  mixtures,  have  been  conducted  over  the  past 
few  years  and  the  results  have  been  recently  reviewed  [1,2].  No 
attempt  will  be  made  to  repeat  those  reviews  here;  however,  it  is 
necessary  to  highlight  certain  studies  of  the  reactivity  of  hydrocar¬ 
bon  mixtures  to  put  the  present  work  in  proper  prospective. 

Leppard  [3]  examined  the  autoignition  of  binary  blends  of 
alkenes/alkanes  in  a  motored  engine  to  investigate  the  chemical 
origins  of  non-linear  blending  octane  numbers.  Stoichiometric  mix¬ 
tures  of  2-butene/n-heptane  and  1,3-butadiene/n-heptane  were  ex¬ 
amined.  The  2-butene/n-heptane  mixture  exhibited  a  significantly 
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higher  octane  quality  than  neat  n-heptane,  but  significantly  lower 
than  neat  2-butene.  Furthermore,  the  intermediates  formed  exclu¬ 
sively  from  n-heptane  occurred  in  the  same  relative  concentrations 
as  observed  for  neat  n-heptane,  however,  this  did  not  apply  to  the 
2-butene  intermediates.  A  similar  phenomena  was  also  observed 
for  the  1,3-butadiene/n-heptane  mixture.  In  the  mixtures  alkenes 
formed  fewer  cyclic  ethers  as  intermediates  than  in  neat  oxidation 
and  instead  formed  aldehydes.  Leppard  suggested  that  the  interac¬ 
tion  between  the  mixture  constituents  occurred  only  through  small 
radicals,  such  as  HO2  and  OH. 

Vanhove  et  al.  [4]  investigated  the  autoignition  of  several  stoi¬ 
chiometric  mixtures  in  a  rapid  compression  machine  below  900  K, 
specifically,  n-heptane/toluene,  iso-octane/toluene,  iso-octane/1 - 
hexene,  1 -hexene/toluene,  and  iso-octane/1 -hexene/toluene  blends. 
During  the  study,  they  measured  ignition  delay  times  and  hydro¬ 
carbon  intermediates  at  selected  residence  times  after  compres¬ 
sion.  Several  significant  differences  between  the  n-heptane/toluene 
and  iso-octane/toluene  blends  were  observed.  The  addition  of 
toluene  to  n-heptane  reduced  the  overall  reactivity  and  the  tem¬ 
perature  dependence  of  the  reactivity.  The  addition  of  the  toluene 
to  iso-octane  altered  the  phenomenology  of  the  autoignition  pro¬ 
cess  by  reducing  the  reactivity  of  the  blend  below  830  K  and 
increasing  the  reactivity  at  higher  temperatures.  The  three  compo¬ 
nent  mixture,  a  volumetric  blend  of  47%  iso-octane,  18%  1 -hexene, 
and  35%  toluene,  produced  cool  flame  behavior.  The  conversion 
of  the  mixture  components  after  the  cool  flame  varied  from  20% 
for  toluene,  30%  for  iso-octane,  and  35%  for  1 -hexene.  The  main 
products  of  the  oxidation  were  consistent  with  the  independent 
reaction  of  the  “pure”  hydrocarbons,  not  from  any  interaction  of 
the  components.  However,  due  to  the  difficulty  in  quantifying  the 
broad  range  of  possible  intermediate  and  final  species  produced 
during  the  oxidation  of  such  a  complex  mixture,  Vanhove  et  al. 
only  accounted  for  approximately  half  of  the  carbon  atoms. 

Khan  [5]  formulated  a  surrogate  mixture,  4.6%  1-pentene,  49.6% 
iso-octane,  14%  n-heptane,  and  31.8%  toluene  by  volume,  to  have 
similar  ratios  of  hydrocarbon  functional  groups,  octane  rating,  and 
low  and  intermediate  temperature  reactivity  as  several  industry 
standard  full  boiling  range  fuels,  specifically  RFA  (Reference  Fuel 
A)  from  the  Auto/Oil  Air  Quality  Improvement  Research  Program, 
indolene  (the  certification  fuel)  and  a  standard  test  fuel  from  Ford. 
The  addition  of  1-pentene  and  toluene  to  a  n-heptane  and  iso¬ 
octane  blend,  reduced  the  overall  reactivity,  shifted  the  start  of 
the  Negative  Temperature  Coefficient  (NTC),  a  region  between  600 
and  800  K  where  reactivity  decreases  as  temperature  increases,  to 
higher  temperatures  in  a  flow  reactor,  and  delayed  the  autoigni¬ 
tion  behavior  in  engines.  A  functional  group  analysis  method  was 
developed  to  use  infrared  absorption  to  provide  semi-quantitative 
measurements  of  aromatics,  alkenes,  aldehydes,  and  alkanes  from 
the  flow  reactor  studies.  Khan  noted  that  the  concentration  of  aro¬ 
matics  was  nearly  constant  and  concluded  that  toluene  did  not 
react  in  the  low  and  intermediate  temperature  regime.  However,  as 
the  infrared  technique  was  tuned  to  aromatic  ring  structure,  these 
results  may  have  simply  indicated  that  the  ring  did  not  open. 

The  present  study  expands  on  this  research  to  elucidate  the  ef¬ 
fects  of  the  alkene  and  aromatic  components  in  full  boiling  range 
fuels  and  to  provide  an  experimental  foundation  for  modeling  sur¬ 
rogate  mixtures.  To  this  end  the  preignition  chemistry  of  the  in¬ 
dividual  components  and  gasoline  surrogate  mixture  developed  by 
Khan,  referred  to  as  GS4c-l,  was  examined  by  measuring  the  re¬ 
acted  gas  composition  at  selected  conditions. 

2.  Experimental  facility 

This  study  utilized  the  Drexel  Pressurized  Flow  Reactor  (PFR) 
Facility  operating  in  the  low  and  intermediate  temperature  regime 
(600-800  K)  at  an  elevated  pressure  (8  atm)  under  dilute  condi¬ 


tions.  The  PFR  facility  was  designed  such  that  chemical  processes 
could  be  examined  in  relative  freedom  from  fluid  mechanics  and 
temperature  gradient  effects.  A  detailed  description  of  the  exper¬ 
imental  facility  and  experimental  techniques  has  been  reported 
earlier  [6-8]  and  only  a  summary  is  presented  here.  The  key  op¬ 
erational  feature  of  the  reactor  is  a  quartz  reactor  tube  within  a 
pressure  vessel.  To  achieve  high  reaction  pressures,  the  inside  of 
the  quartz  reactor  and  the  annulus  between  it  and  the  pressure 
vessel  are  maintained  at  the  same  pressure. 

Nitrogen  (purity  >99.999%)  and  oxygen  (purity  >99.6%)  were 
mixed  into  the  PFR  system  to  create  a  synthetic  “air”  mixture 
free  of  contaminants.  This  “oxidizer”  stream  was  heated  by  in-line 
heaters  to  achieve  the  desired  reaction  temperatures.  To  maximize 
the  temperature  uniformity  along  the  length  of  the  reactor  tube, 
two  techniques  were  utilized.  First,  the  fuel  and  oxidizer  were 
diluted  with  nitrogen  to  limit  the  temperature  rise  due  to  heat  re¬ 
lease.  For  these  experiments,  the  dilution  was  varied  according  to 
the  reactivity  of  the  fuel,  specifically  between  62-85%,  where  the 
dilution  was  defined  as  the  volumetric  flow  of  additional  nitro¬ 
gen  divided  by  the  volumetric  flow  of  the  synthetic  “air”  plus  the 
additional  nitrogen.  Second,  the  walls  of  the  pressure  vessel  were 
heated  by  three  independently  controlled  resistive  bead  heaters. 
Using  this  combination  of  techniques,  the  initial  temperature  pro¬ 
file  was  isothermal  to  within  5°C  without  reaction  and  20  °C  with 
reaction. 

An  independently  controlled  stream  of  heated  nitrogen  (purity 
>99.999%)  was  used  to  vaporize  the  fuel.  The  temperature  of  the 
nitrogen  was  controlled  to  a  temperature  20-50  °C  above  the  fuel 
vaporization  temperature,  as  calculated  from  the  Antoine  equa¬ 
tion  [9].  The  liquid  fuel  was  injected  into  the  centerline  of  the 
nitrogen  stream  by  a  high  pressure  syringe  pump  approximately 
one  meter  from  the  opposed  jet  mixing  nozzle  of  the  reactor.  This 
provided  sufficient  time  to  vaporize  the  fuel  completely  and  to  re¬ 
move  any  inhomogeneities.  The  syringe  pump  permitted  accurate 
(within  0.5%  of  set  point),  repeatable,  and  stable  delivery  of  liquid 
fuels,  enabling  day-to-day  repeatability  of  CO  concentrations  (used 
for  reactivity  mapping)  to  less  than  50  ppm  during  this  study.  To 
create  the  reactive  mixture,  the  oxidizer  stream  and  the  prevap¬ 
orized  fuel/nitrogen  stream  were  rapidly  mixed  in  an  opposed  jet 
annular  mixing  nozzle  at  the  entrance  of  the  quartz  reactor  tube. 

A  water-cooled,  glass-lined  probe  was  used  to  extract  samples 
from  the  centerline  of  the  quartz  reactor  tube  and  to  provide  rapid 
quenching  of  chemical  reactions.  Sample  temperatures  were  mea¬ 
sured  using  a  type-K  thermocouple  integrated  into  the  probe  as¬ 
sembly  and  extending  2.5  cm  past  the  end  of  the  probe  tip.  For 
the  present  experiments,  the  PFR  was  operated  in  a  controlled 
cool  down  (CCD)  mode  to  examine  the  reactivity  over  a  range 
of  temperatures  at  a  fixed  residence  time  and  pressure,  creating 
a  “reactivity  map”  of  the  fuel.  In  this  mode,  the  probe  position 
along  the  centerline  of  the  quartz  tube  is  adjusted  using  a  lin¬ 
ear  positioning  table  in  order  to  maintain  a  fixed  residence  time 
as  the  reactor  slowly  cools  (2-3°C/min).  The  extracted  gas  sam¬ 
ples  were  directed  and  split  among  online  NDIR  (non-dispersive 
infrared)  analyzers  for  carbon  monoxide  (CO)  and  carbon  dioxide 
(CO2)  measurements,  a  total  hydrocarbon  (THC)  analyzer,  and  a 
multiple  loop  sample  (MLS)  storage  system  for  offline  analysis  of 
up  to  14  stored  samples  collected  at  selected  reaction  tempera¬ 
tures.  The  transfer  lines  and  MLS  system  were  heated  to  120  °C  to 
minimize  condensation  of  the  high  boiling  point  hydrocarbons  in 
the  sample,  specifically  higher  molecular  weight  species  and  oxy¬ 
genated  hydrocarbons. 

The  reactivity  of  the  fuel,  including  NTC  behavior,  was  deter¬ 
mined  through  continuous  monitoring  of  CO  formation.  This  mea¬ 
surement  has  been  shown  to  be  a  good  measure  of  low  and  inter¬ 
mediate  temperature  reactivity,  as  CO  does  not  oxidize  significantly 
to  CO2  in  this  temperature  regime  [10]. 
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Table  1 

Summary  of  experimental  conditions. 


Fuel 

Purity 

Equivalence 

ratio 

Dilution 

(%) 

Residence 
time  (ms) 

Fuel  concentration 
(ppm) 

Carbon  atoms 

1-Pentene 

>99% 

0.51 

70 

175  ±  15  ms 

4276 

21380 

n- Heptane 

>99.7% 

0.51 

85 

100  ±7ms 

1459 

10210 

Iso-octane 

>99.9% 

1.0 

62 

250  ±  26  ms 

6346 

50760 

Toluene 

>99.5% 

1.0 

60 

250  ±  26  ms 

6435 

45050 

GS4c-l 

0.73 

65 

225  ±  22  ms 

C5H10:  232 

C7H16:  705 

C8Hi8:  2500 

C7H8:  1602 

Total:  5039 

C5Hi0:  1160 
C7H16:  4935 
C8Hi8:  20000 
C7H8:  11214 

Total:  37309 

600  620  640  660  680  700  720  740  760  780  800 

Temperature  (K) 


Fig.  1.  Reactivity  behavior  of  the  gasoline  surrogate  and  its  components:  (■)  1-Pentene;  (A  )  iso-octane;  (A )  n-heptane;  (•)  surrogate  (GS4c-l). 


The  offline  analysis  of  the  stable  intermediate  species  was  ac¬ 
complished  with  a  gas  chromatograph  (GC)  with  a  flame  ionization 
detector  (FID)  and  a  mass  spectrometer  (MS).  For  chromatographic 
separation,  a  Supelco  Petrocol  DH  column  (100  m,  0.5  pm  film 
thickness,  0.25  mm  OD)  was  utilized.  The  flow  from  the  Petrocol 
DH  column  was  split  between  the  FID  and  MS  detectors  to  per¬ 
mit  simultaneous  species  quantification  and  identification.  Iden¬ 
tification  of  the  unknowns  was  accomplished  by  retention  time 
matching  and/or  by  matching  of  the  measured  MS  spectrum  with 
the  NIST  ’02  database  (145,000+  unique  compounds)  [11].  Con¬ 
centration  of  the  identified  compounds  was  determined  primarily 
using  the  FID  in  the  GC.  Compounds,  for  which  standards  were 
not  available,  were  quantified  employing  the  calibration  curve  of 
a  molecule  with  similar  structure.  Analysis  of  variance  and  regres¬ 
sion  analysis  were  employed  to  determine  the  errors  associated 
with  the  quantification  of  the  intermediate  combustion  species. 
The  results  of  this  analysis  provided  the  uncertainty  in  any  single 
value  calculated  from  the  calibration  curve,  and  these  are  included 
in  the  graphs  for  the  intermediates  measured  with  the  FID. 

3.  Experimental  results 

3.1  Reactivity  mapping 

The  surrogate  and  its  components  were  examined  in  the  low 
and  intermediate  temperature  regime  at  elevated  pressure  (8  atm). 
The  specific  equivalence  ratio  and  dilution  utilized  was  depen¬ 
dent  on  the  reactivity  of  the  fuel,  and  the  specific  values  for  these 
parameters  were  selected  to  limit  heat  release  while  achieving  suf¬ 


ficient  reactivity  for  identification  of  the  major  intermediates.  The 
specific  conditions  for  each  of  the  fuels  are  outlined  in  Table  1. 
Note  that  the  uncertainty  of  residence  time  is  the  uncertainty  in 
the  absolute  sampling  time,  while  the  relative  uncertainty  of  res¬ 
idence  times  within  a  fuel  experiment  is  a  few  milliseconds.  Fur¬ 
thermore,  the  samples  of  the  fuels  were  selected  to  be  as  pure  as 
economically  feasible  to  limit  the  possibility  of  impurities  initiating 
reaction. 

Under  the  conditions  of  the  neat  1-pentene  study,  reactivity 
started  at  640  K  and  CO  concentration  peaked,  indicating  the  start 
of  NTC  behavior,  at  710  K,  Fig.  1.  Neat  n-heptane  started  reacting 
at  630  K  and  also  exhibited  a  strong  NTC  behavior,  observed  to 
begin  at  700  K.  During  the  CCD  experiments  with  2,2,4-trimethyl- 
pentane  (iso-octane),  weak  reactivity  started  at  630  K  and  NTC 
behavior  was  observed  to  begin  at  665  K.  The  40  °C  lower  temper¬ 
ature  for  the  start  of  the  NTC  region  for  iso-octane,  as  compared 
to  n-heptane,  is  intuitively  backward  for  a  fuel  with  a  higher  oc¬ 
tane  rating.  However,  the  presence  of  the  weaker  tertiary  bonded 
hydrogen  in  iso-octane  allowed  the  initiation  to  occur  at  lower 
temperatures.  Neat  toluene  was  not  reactive  at  the  conditions  ex¬ 
amined. 

For  the  surrogate  mixture,  GS4c-l,  reactivity  started  at  630  K 
and  NTC  behavior  was  observed  to  begin  at  693  K.  Although  com¬ 
paring  the  start  of  the  NTC  region  of  the  components  with  the 
surrogate  is  only  strictly  valid  for  identical  conditions,  higher  con¬ 
centrations  should  shift  the  start  of  the  NTC  region  to  higher 
temperatures.  In  these  studies,  the  iso-octane  was  examined  at 
a  higher  concentration  relative  to  the  surrogate  mixture  and  1- 
pentene  was  examined  at  a  lower  concentration.  Considering  the 


552 


D.B.  Lenhert  et  al.  /  Combustion  and  Flame  156  (2009)  549-564 


o 


o 

(/) 

£ 

CD 

L- 

D 

*-» 

5 

CD 

Q. 

E 

CD 


Fig.  2.  Species  profiles  and  temperature  rise  as  a  function  of  reactor  temperature  for  1-pentene  (4276  ppm,  0  =  0.51,  70%  dilution,  175  ms,  and  8  atm):  (■,  — )  CO; 
(A,—  —  )  CO2;  (A, . .)  1-pentene;  .)  AT.  Symbols  are  experimental  results;  lines  are  model  predictions. 


Fig.  3.  Aldehyde  and  ether  species  profiles  as  a  function  of  reactor  temperature  for  1-pentene  (4276  ppm,  </>  =  0.51,  70%  dilution,  175  ms,  and  8  atm):  (■,  — )  Acetaldehyde; 

(A,—  —  )  propanal;  (A, ......)  2-propenal;  (#, -  ■)  butanal;  to-..)  2-pentenal;  (O,  mmmmm)  2-propyl-oxirane.  Symbols  are  experimental  results;  lines  are  model 

predictions. 


concentration  effect,  the  temperature  for  start  of  the  NTC  region 
for  blends  should  be  bound  by  the  start  of  the  NTC  region  for  the 
individual  components  of  that  blend.  A  similar  bounding  was  ob¬ 
served  when  comparing  the  gasoline  primary  reference  fuels  (PRF, 
n-heptane  and  iso-octane)  and  their  blends,  as  the  start  of  the  NTC 
region  for  the  PRF  blends  started  near  687  K  [1].  This  bounding 
suggests  that  the  addition  of  1-pentene  to  the  surrogate  mixture 
could  partially  be  responsible  for  the  shift  of  the  start  of  the  NTC 
region  to  higher  temperatures  as  compared  to  PRF  blends.  How¬ 
ever,  due  to  the  low  concentrations  of  1-pentene  in  the  surrogate 
mixture,  its  more  likely  that  the  combination  of  1-pentene  and 
toluene’s  unreactive  and  radical  scavenging  nature  resulted  in  the 
shift  of  the  start  of  the  NTC  region. 

3.2.  1-Pentene  oxidation  speciation  results 

Samples  were  collected  at  14  temperatures  during  the  1- 
pentene  reactivity  mapping  experiment  to  identify  and  quan¬ 
tify  intermediate  species  in  the  low  and  intermediate  tempera¬ 


ture  regime.  The  identified  intermediates,  signified  by  symbols 
in  Figs.  2-4,  accounted  for  85%  of  the  available  carbon  atoms 
(21380  C  atoms).  Note  that  this  carbon  balance  only  includes  the 
carbon  from  positively  identified  compounds  and  does  not  include 
unidentified  peaks  of  known  molecular  weight.  Including  these 
unidentified  peaks,  which  MS  fragmentation  patterns  suggest  are 
mainly  C5  oxygenates,  would  increase  the  carbon  balance  to  above 
90%. 

Also  shown  in  Figs.  2-4  are  species  profiles  calculated  using  the 
detailed  kinetic  model  of  Touchard  et  al.  [12].  The  mechanism  con¬ 
tained  3385  elementary  reactions  among  837  chemical  species.  In 
this  work,  focusing  on  experimentally  measuring  the  reactivity  and 
reaction  pathways  of  the  selected  hydrocarbons,  no  attempt  was 
made  to  optimize  or  analyze  the  results  of  the  model  simulations. 
The  simulations  are  presented  primarily  to  indicate  the  state  of  our 
chemical  understanding  of  the  reactivity  of  these  hydrocarbons. 

The  PFR  was  modeled  in  Chemkin  3.7.1  as  an  adiabatic  plug 
flow  reactor  40  cm  long  and  2.2  cm  in  diameter  with  no  surface 
chemistry.  The  adiabatic  representation  of  the  reactor  was  chosen 
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Fig.  4.  Alkene  species  profiles  as  a  function  of  reactor  temperature  for  1-pentene  (4276  ppm,  0  =  0.51,  70%  dilution,  175  ms,  and  8  atm):  (M,  ■  )  Ethene;  (A,—  —  ) 

1-propene;  (A, . .)  1, 3-butadiene;  .)  2-pentene;  (Qiw  •  ■ )  cis-  and  trans-l,3-pentadiene.  Symbols  are  experimental  results;  lines  are  model  predictions. 


as  some  temperature  rise  occurred  along  the  length  of  the  reac¬ 
tor.  The  magnitude  of  the  temperature  rise  can  not  be  directly 
measured  during  the  experiments  but  can  be  determined  by  com¬ 
paring  the  time  history  of  the  actual  sample  temperature  to  that  of 
a  “non-reacting”  sample  temperature.  The  “non-reacting”  temper¬ 
ature  history  was  modeled  as  a  second  order  polynomial  and  was 
tuned  to  the  specific  experimental  conditions  by  removing  temper¬ 
atures  where  significant  reactivity  occurred  during  the  experiment. 
Nominally,  this  procedure  permitted  estimation  of  the  temperature 
rise  to  within  a  few  degrees.  Determination  of  this  temperature 
rise  can  assist  in  assuring  that  nearly  isothermal  conditions  are 
achieved  and  aid  in  modeling  the  specific  temperature  conditions 
in  the  flow  reactor.  One  complication  with  modeling  flow  reactor 
data  is  that  there  is  a  finite  physical  mixing  distance  of  the  fuel 
and  oxidizer  streams  and  a  thermal  mixing  distance  due  the  fuel 
and  oxidizer  streams  entering  at  two  different  temperatures,  thus 
resulting  in  a  “poorly”  defined  reaction  time  equals  zero.  As  a  re¬ 
sult,  flow  reactor  model  predictions  are  often  shifted  in  time  to 
match  the  experiments  at  a  specific  condition.  However,  no  time 
shift  was  applied  to  the  simulated  species  profiles  calculated  in 
this  study. 

Fig.  2  presents  the  measured  fuel  mole  fraction  along  with  the 
measurements  of  CO,  CO2,  and  temperature  rise.  Slightly  more  than 
40%  of  the  fuel  was  converted  to  intermediate  species  of  which  CO 
and  CO2  formation  accounted  for  more  than  6%  of  the  total  ini¬ 
tial  carbon  atoms  (21380  C  atoms)  at  the  start  of  the  NTC  region. 
Large  concentrations  of  aldehydes  were  formed  during  the  oxida¬ 
tion  of  1-pentene  with  butanal  and  acetaldehyde  being  the  most 
abundant,  Fig.  3.  While  not  quantified  by  FID,  large  concentrations 
of  formaldehyde  were  detected  with  the  MS.  Approximately  10% 
of  the  total  initial  carbon  atoms  were  converted  to  aldehydes  at 
their  formation  peak,  which  occurred  at  the  same  temperature 
as  the  start  of  the  NTC  region,  710  K.  Inclusion  of  an  estimate 
for  formaldehyde,  using  the  MS  signal,  would  increase  the  over¬ 
all  aldehyde  yield  to  over  16%.  The  temperature  of  the  maximum 
concentration  of  the  individual  species  also  corresponded  to  the 
start  of  the  NTC  region,  except  for  2-pentenal  which  peaked  at 
approximately  690  K.  The  only  positively  identified  ether  was  2- 
propyl-oxirane,  Fig.  3.  Alkene  formation  represented  the  only  other 
functional  group  formed  in  appreciable  amounts,  and  at  their  max¬ 
imum  accounted  for  4%  of  the  total  initial  carbon  atoms.  The 
major  alkene  intermediates  formed,  ethene,  1-propene,  and  1,3- 
pentadiene,  are  presented  in  Fig.  4.  Both  of  the  1,3-pentadiene  iso¬ 


mers  were  quantified,  trans-l,3-pentadiene  was  the  favored  isomer 
producing  twice  the  quantity  of  cis-l,3-pentadiene  over  most  of 
the  temperature  range.  Unlike  the  aldehydes,  the  alkene  formation 
peaked  at  temperatures  higher  than  the  start  of  the  NTC  region. 
The  peak  of  formation  for  the  individual  alkene  species  varied  sig¬ 
nificantly.  Of  the  major  alkenes,  both  ethene  and  1-propene  had 
maximum  formation  around  720  K  and  rapidly  declined  at  higher 
temperatures.  1,3-Pentadiene  peaked  at  even  higher  temperatures, 
approximately  740  K,  and  its  formation  increased  dramatically  at 
the  start  of  the  NTC  region. 

The  results  suggest  that  there  are  two  dominant  pathways  for 
the  decomposition  of  1-pentene,  specifically,  hydrogen  abstraction 
of  the  allylic  hydrogen  and  subsequent  alkane  type  decomposi¬ 
tion  reactions  [13,14]  and  radical  addition  to  the  double  bond 
and  decomposition  through  the  “Waddington”  mechanism  [15,16]. 
As  such,  the  abstraction  of  the  allylic  hydrogen  and  subsequent 
alkane  type  decomposition  was  responsible  for  the  formation  of 
1,3-pentadiene,  ethene,  acetaldehyde,  2-propenal  and  other  inter¬ 
mediates.  If  oxygen  addition  occurred  at  the  allylic  abstraction  site 
followed  by  secondary  hydrogen  abstraction  of  a  primary  hydro¬ 
gen  (to  form  the  alkene  equivalent  of  the  QOOH  radical),  then  the 
preferred  decomposition  pathway  was  through  the  formation  of 
a  lower  molecular  weight  alkene  and  aldehyde,  not  through  the 
formation  of  a  cyclic  ether.  The  prominence  of  this  pathway  ac¬ 
counted  for  the  equal  portions  of  ethene  and  2-propenal  and  low 
production  of  2-ethenyl-oxetane  (less  than  2  ppm)  in  these  exper¬ 
iments.  However,  if  the  secondary  hydrogen  abstraction  removed 
the  secondary  hydrogen  adjacent  to  the  allylic  hydrogen,  then  the 
dominant  decomposition  pathway  was  through  addition  of  a  sec¬ 
ond  oxygen  molecule  to  the  radical  site  and  subsequent  decompo¬ 
sition  to  form  acetaldehyde,  the  third  most  abundant  species.  The 
large  concentrations  of  butanal  and  formaldehyde  suggested  that 
the  Waddington  mechanism  was  the  dominant  pathway  for  the 
decomposition  of  the  fuel.  However,  studies  by  Stark  and  Wadding¬ 
ton  [17]  and  Prahbu  et  al.  [18]  have  questioned  its  importance. 
Prahbu  and  coworkers  observed  an  abundance  of  propanal  and  2- 
methyl-tetrahydrofuran  (THF)  from  the  alkane  type  decomposition 
pathways,  but  did  not  measure  butanal  which  is  expected  from 
the  Waddington  mechanism,  and  thus  argued  that  the  Waddington 
mechanism  was  not  nearly  as  important  as  the  alkane  type  hy¬ 
drogen  abstraction  reactions  for  1-pentene.  However,  propanal  and 
2-methyl-THF  require  the  abstraction  of  a  vinylic  hydrogen,  which 
requires  approximately  6  kcal/mol  more  energy  than  secondary  hy- 
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Fig.  5.  Species  profiles  and  temperature  rise  as  a  function  of  reactor  temperature  for  n-heptane  (1469  ppm,  0  =  0.51,  85%  dilution,  100  ms,  and  8  atm):  C,  -  )  CO; 

(A,—  —  )  CO2;  (A, . .)  n-heptane;  (+, -  ■)  AT.  Symbols  are  experimental  results;  lines  are  model  predictions. 


drogen  bonds.  Prahbu  offered  no  explanation  as  to  why  this  larger 
energy  barrier  alkane  pathway  would  dominate  over  other  lower 
energy  barrier  alkane  pathways.  The  experimental  results  of  the 
present  paper  and  that  of  Minetti  et  al.  [19]  both  show  elevated 
concentrations  of  butanal  with  significantly  lower  concentrations 
of  propanal  and  2-methyl-THF.  This  suggests  that  the  abstraction 
of  the  vinylic  hydrogen  may  not  be  as  important  as  proposed  by 
Prahbu.  Further  discussion  of  these  competitive  reaction  steps  is 
presented  in  the  review  of  Battin-Leclerc  [2]. 

In  general,  the  comparisons  between  the  Touchard  model  and 
experimental  results  were  quite  good.  The  adiabatic  description  for 
the  PFR  agreed  reasonably  well  with  the  experimental  temperature 
rise,  albeit  slightly  overpredicting  the  temperature  rise  between 
680  and  750  K,  Fig.  2.  However,  the  results  do  clearly  demon¬ 
strate  why  the  adiabatic  model  was  chosen  for  this  experiment 
over  an  isothermal  model.  The  relatively  minor  overprediction  of 
the  temperature  rise  in  the  adiabatic  representation  of  the  reactor 
can  partially  explain  the  observed  differences  in  fuel  decomposi¬ 
tion  between  the  experiment  and  modeling.  However,  where  the 
temperatures  matched  (<680  K),  sizable  differences  in  intermedi¬ 
ates  were  observed.  Specifically,  the  formation  of  CO,  1-propene, 
and  butanal  were  under  predicted  and  1,3-butadiene  was  over 
predicted.  The  model  also  does  not  predict  the  rapid  increase  in 
1,3-pentadiene  formation  as  the  temperature  increases  beyond  the 
start  of  the  NTC  region.  As  noted  previously,  the  formation  of  1,3- 
pentadiene  occurs  through  alkane  type  decomposition  pathways. 
For  n-heptane,  a  similar  rapid  increase  was  observed  for  the  C7 
alkenes  and  ethers  and  attributed  to  the  backward  shift  for  the 
addition  reaction  of  O2  to  heptyl  and  heptylhydroperoxy  radicals 
(see  below),  suggesting  that  improvements  could  be  made  to  the 
alkane  type  decomposition  pathways  in  the  1-pentene  mechanism. 
Furthermore,  the  mechanism  tended  to  underestimate  the  forma¬ 
tion  of  intermediates  through  the  “Waddington”  reaction  pathway. 
Increasing  the  importance  of  the  “Waddington”  mechanism  ver¬ 
sus  hydrogen  abstraction  of  the  allylic  hydrogen  may  improve  the 
model  predictions  presented  here  and  in  Minetti  et  al.  [19]. 

As  the  alkene  molecules  become  larger  and  there  is  larger  sep¬ 
aration  between  the  double  bond  and  the  terminal  methyl  group, 
it  is  reasonable  to  assume  radical  addition  to  the  double  bond 
would  become  less  important.  It  is  important  to  establish  the  car¬ 
bon  chain  length  that  produces  the  switch  between  the  radical 
addition  pathway  and  the  alkane  hydrogen  abstraction  pathways  in 


order  to  allow  the  simplification  of  the  alkene  decomposition  path¬ 
ways  in  the  surrogate  fuel  mechanisms.  The  results  presented  here 
provide  important  experimental  evidence  as  to  the  importance  of 
the  “Waddington”  mechanism  for  alkenes  in  surrogate  mixtures. 

3.3.  n-Heptane  oxidation  speciation  results 

The  major  intermediates  were  measured  in  samples  extracted 
at  14  temperatures  during  the  n-heptane  reactivity  mapping  ex¬ 
periment.  The  identified  intermediates,  signified  by  symbols  in 
Figs.  5-8,  accounted  for  60%  of  the  available  carbon  atoms  (10210  C 
atoms).  MS  fragmentation  patterns  indicate  that  the  unidentified 
compounds  were  likely  oxygenated  intermediates,  but  they  were 
not  positively  identified.  Positive  identification  of  these  larger  oxy¬ 
genated  compounds  are  very  difficult,  as  they  are  not  available  for 
purchase  nor  are  they  included  in  the  NIST  MS  database.  Using  the 
known  molecular  weights  of  the  next  5  major  intermediates  would 
increase  the  carbon  balance  to  over  80%.  As  with  1-pentene  sim¬ 
ulations  using  an  existing  detailed  kinetic  model  are  also  plotted. 
The  basic  mechanism  and  thermochemistry  by  Curran  et  al.  [13], 
2450  elementary  reactions  among  550  chemical  species,  includ¬ 
ing  updated  chemical  rate  and  thermodynamic  parameters  [20], 
was  utilized  for  this  comparison.  Again,  the  PFR  was  modeled  in 
Chemkin  as  an  adiabatic  plug  flow  reactor  40  cm  long  with  a 
2.2  cm  diameter  with  no  time  shift  or  surface  chemistry. 

Fig.  5  illustrates  the  measured  n-heptane  disappearance  along 
with  the  measurements  of  CO,  CO2,  and  temperature  rise  as  a  func¬ 
tion  of  sample  temperature.  Over  80%  of  the  initial  fuel  present 
was  converted  to  intermediate  species  near  the  start  of  the  NTC 
region,  of  which  CO  and  CO2  formation  accounted  for  more  than 
13%  of  the  total  initial  carbon  atoms  (10210  C  atoms).  Ethers  were 
the  most  abundant  functional  group  measured  over  much  of  the 
temperature  range,  Fig.  6.  At  temperatures  lower  than  the  start 
of  the  NTC  region  ether  formation  was  roughly  constant  and  ac¬ 
counted  for  8%  of  the  total  initial  carbon  atoms.  After  the  start  of 
the  NTC  region,  ether  formation  increased  considerably  to  its  peak 
of  over  16%  near  750  K.  However,  there  were  a  significant  number 
of  ethers  predicted  by  the  mechanism  that  were  not  identified.  The 
major  ether  identified  was  2-methyl-5-ethyl-THF,  Fig.  6.  Both  of  the 
2-methyl-5-ethyl-THF  isomers  were  quantified,  trans-2-methyl-5- 
ethyl-THF  was  the  favored  isomer  producing  1.6  times  the  quantity 
of  cis-isomer  prior  to  the  start  of  the  NTC  region  and  1.3  times 
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Fig.  6.  Alkene  and  ether  species  profiles  as  a  function  of  reactor  temperature  for  n-heptane  (1469  ppm,  0  =  0.51,  85%  dilution,  100  ms,  and  8  atm):  (M,  )  Ethene; 

(A ,—  —  )  1-propene;  (A, . ,)  2-propyl-THF;  (#,^—  .)  2-methyl-5-ethyl-THF;  (Q^i  1 1 )  2-methyl-4-propyl-oxetane.  Symbols  are  experimental  results;  lines  are  model 

predictions. 


Fig.  7.  Alkene  profiles  as  a  function  of  reactor  temperature  for  n-heptane  (1469  ppm,  0  =  0.51,  85%  dilution,  100  ms,  and  8  atm):  (■,  ■  )  1-Butene;  (A  ,  —  —  )  1-pentene; 

(A, ......)  1 -hexene;  (•,—  .)  1-heptene;  (D,™  <  i )  cis-  and  trans-2-heptene;  (O, wihmhwwii)  trans-3-heptene.  Symbols  are  experimental  results;  lines  are  model  predictions. 


in  the  NTC  region.  Alkenes  (Figs.  6  and  7)  were  the  second  most 
abundant  functional  group  over  much  of  the  experimental  condi¬ 
tions,  representing  nearly  12%  of  the  total  initial  n-heptane  carbon 
atoms  at  its  maximum  with  ethene  and  1-propene  being  the  most 
abundant.  The  alkene  formation  peaked  at  temperatures  above 
the  start  of  the  NTC  region,  roughly  40  °C  higher,  similar  to  the 
ethers  and  1-pentene  results.  The  temperature  of  peak  formation 
for  the  individual  alkene  species  did  not  vary  as  much  as  in  the 

1- pentene  oxidation.  Most  of  the  major  alkenes  peaked  between 
740  and  750  K,  the  most  notable  exception  was  ethene  which 
peaked  around  720  K.  Almost  all  of  conjugate  alkenes  were  identi¬ 
fied,  the  only  exception  was  cis-3-heptene,  because  it  eluted  with 
n-heptane  during  GC  analysis.  Based  on  the  measurements  of  the 
individual  isomers  of  2-heptene,  trans-2-heptene  was  the  favored 
isomer  resulting  in  a  2:1  ratio  over  cis-2-heptene.  Comparing  the 

2- heptene  isomers  to  1-heptene  showed  that  the  2-heptene  iso¬ 
mers  were  favored  and  the  degree  of  preference  increased  steadily 
from  1.2  to  1.6  as  the  temperature  increased  to  750  K.  Furthermore, 
comparing  trans-2-heptene  with  trans-3-heptene  showed  that  the 


3-heptene  isomer  was  favored  slightly  at  a  1.3:1  ratio.  This  sug¬ 
gests  that  overall  the  3-heptene  isomer  was  the  preferred  conju¬ 
gate  alkene  with  2-heptene  being  next  most  preferred.  However, 
this  was  expected  as  the  formation  of  1-heptene  requires  the  ab¬ 
straction  of  a  stronger  primary  hydrogen  bond.  A  significant  fea¬ 
ture  of  the  conjugate  alkenes,  and  that  of  the  C7  ethers,  was  that 
their  concentration  rose  significantly  as  the  temperature  increased 
beyond  the  start  of  the  NTC  region.  This  sudden  increase  coin¬ 
cided  with  a  decrease  in  aldehyde  concentration.  These  shifts  were 
due  to  the  backward  shift  of  O2  addition  to  heptyl  (C7H15)  and 
heptylhydroperoxy  (C7H14OOH)  radicals  [21,22].  As  the  equilibrium 
shifts  backward  at  the  start  of  the  NTC  region,  a  larger  fraction  of 
C7H15  and  C7H14OOH  radicals  were  available  to  undergo  alterna¬ 
tive  reactions.  Hence,  more  conjugate  alkenes  are  formed  from  the 
decomposition  of  C7H15  and  more  ethers  are  formed  from  the  de¬ 
composition  of  C7H14OOH. 

Aldehyde  formation  ranged  from  8  to  9%  of  the  total  initial 
carbon  atoms  for  much  of  the  temperature  range  examined  with 
acetaldehyde  and  propanal  reaching  their  peak  at  approximately 
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Fig.  8.  Aldehyde  profiles  as  a  function  of  reactor  temperature  for  n-heptane  (1469  ppm,  0  =  0.51,  85%  dilution,  100  ms,  and  8  atm):  (M,  ■  )  Acetaldehyde:  (A,—  —  ) 

propanal;  (A, . .)  2-propenal;  (•,—  .)  butanal;  i )  pentanal.  Symbols  are  experimental  results;  lines  are  model  predictions. 
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Fig.  9.  Structure  and  H-abstraction  sites  of  n-heptane  and  iso-octane. 


680  K,  Fig.  8.  However,  inclusion  of  a  formaldehyde  production  es¬ 
timate  from  the  MS  would  increase  the  overall  aldehyde  yield  to 
approximately  14%  and  the  formation  would  peak  at  a  tempera¬ 
ture  near  the  start  of  the  NTC  region.  Unlike  the  aldehyde  forma¬ 
tion  during  1-pentene  oxidation,  the  individual  aldehyde  species 
formed  during  n-heptane  oxidation  were  not  as  sensitive  to  tem¬ 
perature,  as  the  profiles  essentially  rose  to  a  plateau  level  over  the 
range  of  experimental  conditions.  From  a  mechanistic  perspective, 
two  major  pathways  have  been  used  to  describe  the  formation  of 
the  aldehydes,  specifically  through  the  decomposition  of  the  QOOH 
and  OOQOOH  radicals.  The  decomposition  of  the  C7H14OOH  radi¬ 
cal  would  result  in  equal  concentrations  of  lower  molecular  weight 
alkene  and  aldehyde  pairs,  for  example  1 -butene  and  propanal  or 
1-pentene  and  acetaldehyde.  However,  the  concentrations  of  these 
aldehyde  and  alkene  pairs  were  not  measured  in  equal  proportions. 
This  suggests  that  the  branching  pathway  decomposing  the  dihy- 
droperoxyheptyl  radical  (OOC7H14OOH)  was  the  primary  source  of 
the  aldehydes.  Therefore,  the  two  most  abundant  aldehydes,  ac¬ 
etaldehyde  and  propanal,  resulted  from  the  decomposition  of  a 
bC7Hi300-b00H  radical  and  a  CC7H13OO-COOH  radical,  respec¬ 
tively,  where  the  letters  “b”  and  “c”  indicate  a  distinct  site  for 
hydrogen  abstraction,  Fig.  9.  The  only  exception  for  this  dominant 
pathway  for  aldehyde  formation  is  butanal.  The  equal  proportions 
of  butanal  and  1-propene,  the  6th  and  7th  most  abundant  species, 
indicate  their  formation  was  a  result  of  the  decomposition  of  the 
dC7Hi4OOH-b  radical. 

In  general,  the  model  predictions  of  n-heptane  reaction  were 
generally  good,  however,  at  temperatures  above  670  K,  the  mech¬ 
anism  overpredicted  the  extent  of  fuel  conversion.  As  with  1- 
pentene,  the  temperature  rise  with  the  adiabatic  representation  of 


the  reactor  was  higher  than  the  experiments.  After  the  start  of 
the  NTC  region,  the  slope  of  the  experimental  and  model  tem¬ 
perature  rise  was  nearly  identical,  however,  the  model  predic¬ 
tions  for  fuel  concentrations  diverged  significantly  from  the  exper¬ 
iments.  Specifically,  the  formation  of  the  conjugate  alkenes,  Fig.  7, 
was  underpredicted  by  approximately  a  factor  of  two,  despite  the 
overprediction  of  fuel  disappearance.  The  most  significant  differ¬ 
ences  were  observed  in  the  lower  molecular  weight  aldehydes  and 
alkenes,  namely  ethene,  acetaldehyde,  and  propanal.  Nevertheless, 
the  general  trends  and  temperature  dependence  of  the  intermedi¬ 
ate  species  were  captured  reasonably  well.  Perhaps  most  notable  is 
the  proper  calculation  of  the  trends  for  the  ethers,  Fig.  6,  as  there 
is  a  lack  of  data  available  to  validate  submechanisms  describing 
their  decomposition. 

3.4.  Iso-octane  oxidation  speciation  results 

Samples  were  only  extracted  and  analyzed  at  13  temperatures 
during  the  iso-octane  (2,2,4-trimethyl-pentane)  experiment.  The 
identified  intermediates,  signified  by  symbols  in  Figs.  10-13,  ac¬ 
counted  for  100  ±  3%  of  the  available  carbon  atoms  (50760  C 
atoms).  Several  lower  concentration  intermediates,  likely  oxy¬ 
genated  intermediates,  were  not  identified.  If  included  in  the  car¬ 
bon  balance,  they  would  not  have  had  a  significant  impact  on 
the  overall  carbon  balance.  As  with  n-heptane,  the  intermediate 
species  profiles  were  compared  to  predictions  using  an  updated 
version  of  a  detailed  iso-octane  mechanism  published  by  Curran  et 
al.  [20,23]  containing  3600  elementary  reactions  among  860  chem¬ 
ical  species.  The  higher  initial  iso-octane  concentration  along  with 
the  overprediction  of  the  iso-octane  reactivity  from  the  model,  re- 
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Fig.  10.  Species  profiles  and  temperature  rise  as  a  function  of  reactor  temperature  for  iso-octane  (6346  ppm,  0  =  1.0,  62%  dilution,  250  ms,  and  8  atm):  (■,  — )  CO; 
(A,—  -  )  C02;  (♦, . >)  iso-octane;  (+,—  •)  AT.  Symbols  are  experimental  results;  lines  are  model  predictions. 


Fig.  11.  Oxygenate  species  profiles  as  a  function  of  reactor  temperature  for  iso-octane  (6346  ppm,  0  =  1.0,  62%  dilution,  250  ms,  and  8  atm):  (■,  — )  Acetaldehyde; 

(A,—  —  )  2-methyl-propanal;  (A, . .)  2,2-dimethyl-propanal;  (•,—  -)  2-propanone;  (dH,i™  1 i)  2-iso-propyl-3,3-dimethyl-oxetane;  (O, 2,2,4,4-tetramethyl- 

THF.  Symbols  are  experimental  results;  lines  are  model  predictions. 


suited  in  excessive  heat  release  and  temperature  rise  predictions 
using  the  adiabatic  representation  of  the  reactor.  As  a  result,  the 
isothermal  representation  was  used  as  it  provided  a  more  realistic 
depiction  of  the  experimental  temperature  profiles. 

Fig.  10  shows  the  iso-octane  mole  fractions  along  with  the 
measurements  of  CO,  CO2,  and  temperature  rise  as  a  function  of 
sample  temperature.  Only  20%  of  the  fuel  was  converted  to  in¬ 
termediate  species  at  the  start  of  the  NTC  region.  CO  and  CO2 
formation  accounted  for  less  than  1%  of  the  total  initial  carbon 
atoms  (50760  C  atoms)  at  the  start  of  the  NTC  region.  The  major 
class  of  intermediates  formed  before  the  start  of  the  NTC  region 
was  ethers  which  account  for  over  11%  of  the  total  initial  car¬ 
bon  atoms.  The  ether  2,2,4,4-tetramethyl-THF  accounted  for  over 
65%  of  the  identified  ethers,  Fig.  11.  A  similar  preferential  forma¬ 
tion  of  a  single  ether  was  also  observed  with  n-heptane.  Only  a 
few  aldehydes  and  ketones  were  identified  with  2-propanone  and 
acetaldehyde  being  the  most  abundant,  Fig.  11.  Aldehydes  and  ke¬ 
tones  accounted  for  less  than  2%  and  1%  of  the  total  initial  car¬ 
bon  atoms  at  the  maximum,  respectively.  Inclusion  of  a  formalde¬ 
hyde  production  estimate  would  only  increase  the  overall  aldehyde 
yield  to  slightly  more  than  2%.  Alkenes  (Figs.  12  and  13)  repre¬ 


sented  the  second  most  abundant  functional  group  to  be  formed 
over  this  temperature  range  and  represented  4%  of  the  total  ini¬ 
tial  carbon  atoms  at  its  maximum  with  2-methyl-l-propene  and 
2,4,4-trimethyl-l-pentene  being  the  most  abundant.  The  peak  of 
the  individual  alkene  species  varied  slightly,  but  most  of  the  ma¬ 
jor  alkenes  peaked  at  roughly  670  K,  the  most  notable  exceptions 
were  the  methyl-1 -pentenes  that  peaked  around  660  K.  Both  of 
the  4,4-dimethyl-2-pentene  isomers  were  observed  with  the  cis- 
isomer  formed  at  nearly  3.5  times  the  level  of  the  trans-isomer. 

Although  n-heptane  and  iso-octane  are  both  alkanes,  they  show 
significant  differences  in  the  preferred  decomposition  pathways. 
First,  iso-octane  produced  significantly  lower  concentrations  of 
aldehydes  as  compared  to  n-heptane.  The  dominant  pathway  for 
aldehyde  formation  with  n-heptane  required  two  internal  hydro¬ 
gen  atom  isomerizations.  These  isomerizations  reactions  are  signif¬ 
icantly  slower  in  iso-octane  due  to  its  branched  nature,  resulting  in 
the  significantly  lower  concentrations  of  aldehydes.  In  n-heptane, 
the  other  dominate  pathway  for  aldehyde  formation  occurred  from 
the  decomposition  of  the  QOOH  radical.  This  pathway  tends  to  be 
more  important  in  linear  alkanes  when  the  two  abstraction  sites 
are  2  carbon  atoms  away  from  each  other,  as  indicated  by  the  sig- 
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Fig.  12.  Major  alkene  species  profiles  as  a  function  of  reactor  temperature  for  iso-octane  (6346  ppm,  0  =  1.0,  62%  dilution,  250  ms,  and  8  atm):  (M,  ■  )  2,2,4-Trimethyl- 

1-pentene;  (A ,  —  —  )  2,2,4-trimethyl-2-pentene;  (♦, . •)  2-methyl-  1-propene.  Symbols  are  experimental  results;  lines  are  model  predictions. 


Fig.  13.  Alkene  species  profiles  as  a  function  of  reactor  temperature  for  iso-octane  (6346  ppm,  0  =  1.0,  62%  dilution,  250  ms,  and  8  atm):  (■,  — )  2,4-Dimethyl- 1-pentene; 

(A ,—  —  )  2,4-dimethyl-2-pentene;  (♦, . .)  4,4-dimethyl- 1-pentene;  (#,^—  .)  4,4-dimethyl-2-pentene;  (ED,*  i  ■ )  1-propene.  Symbols  are  experimental  results;  lines 

are  model  predictions. 


nificant  formation  of  butanal  resulting  from  the  decomposition  of 
the  dCyH^OOH-b  radical.  This  would  suggest  that  the  decomposi¬ 
tion  of  the  aCsHi600H-b  radical  would  tend  to  form  an  aldehyde. 
However,  due  to  the  branched  structure  of  iso-octane,  the  pre¬ 
ferred  pathway  was  through  the  formation  of  an  ether,  specifically 
2-isopropyl-3,3-dimethyl-oxetane.  The  consequence  of  these  differ¬ 
ences  resulted  in  the  reduced  importance  of  aldehyde  formation 
from  iso-octane  and  an  associated  reduction  in  CO  formation.  Sec¬ 
ond,  the  relative  importance  of  the  individual  conjugate  alkenes 
for  iso-octane  were  not  ranked  according  to  their  activation  en¬ 
ergies.  The  formation  of  the  conjugate  alkenes,  2, 4, 4-trimethyl- 1- 
pentene  and  2,4,4-trimethyl-2-pentene,  occurred  from  the  internal 
isomerization  of  a  cCsHnOO  radical.  Either  the  “b”  or  “d”  hy¬ 
drogen  must  be  abstracted  to  form  the  conjugate  alkenes,  since 
the  “b”  hydrogen  is  a  secondary  bonded  hydrogen  one  would  ex¬ 
pect  2,4,4-trimethyl-2-pentene  to  be  the  favored  species.  However, 
three  times  the  number  of  “d”  hydrogen  atoms  compared  to  “b” 
hydrogen  atoms  resulted  in  the  higher  observed  concentrations  of 
2,4,4-trimethyl-l-pentene  over  that  of  2,4,4-trimethyl-2-pentene. 


In  general,  the  model  predictions  were  very  poor,  significantly 
over  estimating  the  decomposition  of  the  fuel  and  formation  of 
many  of  the  intermediates.  Specifically,  the  model  overpredicted 
the  formation  of  several  lower  molecular  weight  alkenes  and  oxy¬ 
genates,  namely  2-methyl- 1-propene,  1-propene,  and  2-propanone. 
Due  to  these  significant  discrepancies  between  model  and  ex¬ 
periment,  the  mechanism  was  examined  to  isolate  pathways  and 
intermediates  which  could,  in  part,  account  for  the  observed  dif¬ 
ferences.  The  results  suggest  that  the  rate  of  secondary  oxygen  ad¬ 
dition  to  the  CgHieOOH  radical  may  be  overestimated.  As  a  result, 
radicals  such  as  aCgHieOOH-c  and  aCgHieOOH-b  decompose  to  2- 
propanone  and  2-methyl-propanal  rather  than  2,2,4,4-tetramethyl- 
THF  and  2-isopropyl-3,3-dimethyl-oxetane,  respectively.  This  be¬ 
havior  would  account  for  the  overestimation  of  2-propanone  and 
2-methyl-propanal  and  underestimation  of  2,2,4,4-tetramethyl-THF 
and  2-isopropyl-3,3-dimethyl-oxetane.  Consequently,  the  additional 
dCgHigOOH-c  and  bCgHigOOH-c  radicals  would  result  in  an  in¬ 
crease  in  2,2,4-trimethyl-l-pentene  and  2,2,4-trimethyl-2-pentene, 
which  were  also  underestimated  by  the  mechanism. 
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Fig.  14.  Species  profiles  and  temperature  rise  versus  temperature  for  GS4c-l  (5039  ppm,  (j)  —  0.73,  65%  dilution,  225  ms,  and  8  atm):  (M,  ■■  )  CO;  (A,—  —  )  CO2; 

(A, . .)  1-pentene;  (•,—  .)  n-heptane;  (Q™  1 1)  iso-octane;  (O)  toluene,  and  (“I")  AT.  Symbols  are  experimental  results;  lines  are  model  predictions. 


3.5.  Surrogate  mixture  oxidation  speciation  results 

Reaction  gas  samples  were  collected  at  14  temperatures  during 
the  experiments  with  GS4c-l.  The  identified  intermediates,  signi¬ 
fied  by  symbols  in  Figs.  14,  16-21,  accounted  for  better  than  90% 
of  the  available  carbon  atoms  (37309  C  atoms).  This  carbon  balance 
only  includes  the  carbon  from  positively  identified  compounds  and 
does  not  include  unidentified  peaks  of  known  molecular  weight. 
Including  these  unidentified  peaks  was  estimated  to  account  for  an 
additional  5-7%  of  the  carbon  atoms.  A  detailed  chemical  model 
which  is  able  to  model  the  surrogate  mixture  over  the  condi¬ 
tions  of  this  experiment  was  not  available  in  the  open  literature. 
A  particular  problem  is  that  there  is  a  complete  lack  of  a  toluene 
mechanism  for  lower  temperatures  [2].  Therefore,  the  intermedi¬ 
ate  species  profiles  were  compared  to  predictions  using  a  detailed 
chemical  model  for  n-heptane/iso-octane  oxidation  published  by 
Curran  et  al.  [20,24]  containing  4200  elementary  reactions  among 
1050  chemical  species  and  which  includes  a  basic  1-pentene  sub¬ 
mechanism.  This  mechanism  was  recently  used  to  model  homo¬ 
geneous  charge  compression  ignition  (HCCI)  engines  with  reason¬ 
able  success  [25].  The  initial  conditions  for  the  surrogate  mixture 
were  modified  to  include  only  the  mole  fractions  of  n-heptane, 
iso-octane,  and  1-pentene,  with  nitrogen  substituted  for  the  mole 
fraction  of  toluene.  This  substitution  allowed  the  isolation  of  the 
effects  of  toluene  addition  on  the  surrogate  mixture,  discussed  be¬ 
low.  As  with  iso-octane,  the  isothermal  representation  was  used  as 
it  provided  a  more  realistic  depiction  of  the  experimental  temper¬ 
ature  profiles  as  the  kinetic  model  significantly  overpredicted  fuel 
consumption  and  energy  release. 

Fig.  14  presents  the  measured  fuel  mole  fractions  and  measure¬ 
ments  of  CO,  CO2,  and  temperature  rise,  along  with  the  associated 
model  predictions.  The  fractional  consumption  of  the  individual 
constituents  in  the  surrogate  varied  significantly,  Fig.  15,  and  is 
related  to  their  respective  research  octane  numbers  (RON).  Over 
60%  of  the  n-heptane  (RON  =  0)  was  converted  to  intermediates, 
while  only  45%  of  the  iso-octane  (RON  =  100)  was  converted. 
The  peak  of  fuel  consumption  occurred  at  the  start  of  the  NTC 
region  for  both  of  these  fuels  and  coincided  with  the  overall  max¬ 
imum  in  fuel  consumption.  The  measurement  of  1-pentene  (RON 
=  91)  included  the  1-pentene  produced  from  n-heptane  oxidation. 
To  correct  for  the  effects  of  this  production,  an  estimate  of  the 
amount  of  1-pentene  produced  from  n-heptane  was  based  upon 
the  ratio  of  1-pentene  to  n-heptane  for  the  neat  fuel  study.  Scaling 


this  ratio  to  the  surrogate  mixture  concentrations  suggested  that  a 
maximum  of  12  ppm,  or  approximately  5%  of  the  initial  1-pentene, 
would  be  produced  from  n-heptane  reaction  in  the  surrogate  mix¬ 
ture.  Furthermore,  the  maximum  formation  of  1-pentene  occurred 
at  a  temperature  higher  than  the  start  of  the  NTC  region  dur¬ 
ing  the  neat  n-heptane  study.  This  explains  why  the  maximum 
percent  consumption  of  1-pentene  shown  in  Fig.  15  is  at  a  tem¬ 
perature  lower  than  the  start  of  the  NTC  region,  680  K.  Toluene 
(RON  =  120)  measurements  exhibited  an  anomalous  behavior  dur¬ 
ing  the  CCD  experiments  and  during  separate  Constant  Inlet  Tem¬ 
perature  (CIT)  experiments  (not  shown  here).  During  the  fuel  cal¬ 
ibration,  toluene  would  be  measured  in  the  proper  percentages 
based  upon  the  surrogate  composition  and  reaction  conditions,  but 
when  samples  were  quantified  under  an  oxidizing  environment 
approximately  10%  more  toluene  was  present  than  expected.  Al¬ 
though  some  of  the  intermediates  from  n-heptane  and  iso-octane 
would  coelute  with  toluene  at  conditions  of  high  n-heptane  and 
iso-octane  reactivity,  these  intermediates  are  not  significant  under 
low  reactivity  conditions,  specifically,  the  600  K  and  780  K  sample 
points,  where  the  difference  was  also  observed.  It  is  unclear  what 
caused  the  variation  in  observed  concentrations.  Despite  this  dis¬ 
crepancy,  the  difference  still  shows  that  over  14%  of  the  toluene 
was  consumed  and  the  peak  reactivity  occurred  at  the  start  of  the 
NTC  region. 

Slightly  more  than  35%  of  the  surrogate  fuel  was  converted  to 
intermediate  species  of  which  CO  and  CO2  formation  accounted 
for  approximately  3%  of  the  total  initial  carbon  atoms  (37309  C 
atoms)  at  the  start  of  the  NTC  region.  Ethers  represented  the  most 
abundant  functional  group  formed  over  the  temperature  range  in¬ 
vestigated  and  represented  a  maximum  of  approximately  10%  of 
the  total  initial  carbon  atoms,  Fig.  16.  The  ethers  attributable  to 
n-heptane  decomposition  were  2-methyl-5-ethyl-tetrahydrofuran 
(THF)  and  2-propyl-THF,  and  those  attributable  to  iso-octane  de¬ 
composition  were  2,2,4,4-tetramethyl-THF  and  2-isopropyl-3,3- 
dimethyl-oxetane.  The  major  ether  formed  by  1-pentene,  2-propyl- 
oxirane,  was  not  quantified  as  it  coeluted  with  iso-octane.  The 
peak  in  the  rate  of  formation  for  the  2-methyl-5-ethyl-THF  iso¬ 
mers  occurred  at  720  K,  approximately  30  °C  after  the  start  of  the 
NTC  region.  The  maximum  concentration  of  both  of  the  Cs  ethers 
coincided  with  the  start  of  the  NTC  region. 

Aldehyde  formation  peaked  near  5%  of  the  total  initial  carbon 
atoms,  the  major  intermediates  are  shown  in  Figs.  17  and  18.  While 
not  quantified  by  FID,  large  concentrations  of  formaldehyde  were 
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Fig.  15.  Overall  fuel  conversion  and  conversion  of  the  individual  components  for  GS4c-l:  (■)  1-Pentene;  (A)  n-heptane;  (A)  iso-octane;  (•)  toluene;  (E3)  overall  fuel 
conversion;  (O)  carbon  balance. 


Fig.  16.  Ether  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  </>  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  C,  )  2-Methyl-5-ethyl-THF; 
(A,—  —  )  2,2,4,4-tetramethyl-THF;  (A, . .)  2-isopropyl-3,3-dimethyl-oxetane;  and  (•,—  .)  2-propyl-THF.  Symbols  are  experimental  results;  lines  are  model  predic¬ 

tions. 


detected  with  the  MS.  Using  the  MS  signal  and  associated  cali¬ 
bration  curve  for  propanal;  the  concentration  of  formaldehyde  can 
be  estimated.  Inclusion  of  this  estimate  would  increase  the  overall 
aldehyde  yield  to  approximately  6.5%.  The  aldehydes  attributable 
to  the  reaction  of  n-heptane,  iso-octane,  and  1-pentene  were  the 
most  abundant,  specifically,  formaldehyde  and  acetaldehyde.  More 
interesting  are  measurements  of  aldehydes  specific  to  a  particular 
fuel  component,  notably,  iso-octane  reaction  produced  2-methyl- 
propanal  and  2,2-dimethyl-propanal  and  n-heptane  produced  pen- 
tanal.  The  other  linear  aldehydes  are  possibly  produced  from  both 
n-heptane  and  1-pentene  consumption. 

Alkene  formation  represented  the  third  most  abundant  func¬ 
tional  group  over  much  of  the  temperature  regime  accounting  for 
more  than  4%  of  the  total  initial  carbon  atoms  at  its  maximum, 
Figs.  19  and  20.  The  alkene  formation  peaked  at  a  temperature 
10  °C  higher  than  the  start  of  the  NTC  region.  The  peak  of  forma¬ 
tion  of  the  individual  alkene  species  varied,  but  most  of  the  major 
alkenes  peaked  at  or  above  the  start  of  the  NTC  region.  Similar  to 
neat  iso-octane,  2,4,4-trimethyl-l-pentene  in  the  surrogate  mixture 


peaked  prior  to  the  start  of  the  NTC  region,  680  K.  Significant  con¬ 
centrations  of  alkenes  specific  to  a  particular  fuel  component  were 
identified,  particularly,  iso-octane  produced  2-methyl-l-propene, 
2,4,4-trimethyl-l-pentene,  and  2,4,4-trimethyl-2-pentene  and  n- 
heptane  produced  cis-  and  trans-2-heptene,  and  trans-3-heptene. 
Intermediates  produced  by  1-pentene,  notably  1,3-butadiene  and 
1,3-pentadiene,  were  only  identified  at  trace  levels.  Vanhove  and 
coworkers  [4]  noted  that  1 -hexene  produced  a  relatively  larger 
concentration  of  hexadiene  for  their  surrogate  mixture.  However, 
due  to  the  significantly  lower  concentration  of  alkene  used  for  the 
surrogate  in  this  study,  this  phenomena  was  not  observed  for  1- 
pentene. 

During  the  oxidation  of  the  surrogate,  three  major  species  were 
quantified  that  could  only  be  produced  from  toluene,  specifically, 
benzaldehyde,  benzene,  and  phenol,  Fig.  21.  Benzaldehyde  repre¬ 
sented  the  most  abundant  ring  species  accounting  for  over  7%  of 
the  total  initial  carbon  atoms  in  toluene  (11214  C  atoms).  Phe¬ 
nol  and  benzene  were  produced  in  nearly  identical  concentra¬ 
tions  and  accounted  for  approximately  1%  of  the  total  initial  car- 
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Fig.  17.  Aldehydes  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  0  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  (M,  -  )  Acetaldehyde: 

(A,—  —  )  propanal;  (A, . .)  2-propenal;  .)  2-methyl-propanal;  and  (EH,™  ■  i)  2-methyl-2-propenal.  Symbols  are  experimental  results;  lines  are  model  pre¬ 

dictions. 


Fig.  18.  Aldehydes  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  0  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  (■,  — )  Butanal;  (A,—  —  ) 
pentanal;  and  (A, . .)  2,2-dimethyl-propanal.  Symbols  are  experimental  results;  lines  are  model  predictions. 


bon  atoms  in  toluene.  Due  to  the  complexity  of  the  gas  mixture 
and  the  possibility  of  coelution  of  toluene  reaction  intermedi¬ 
ates,  the  chromatograms  were  scrutinized  for  additional  aromat¬ 
ics  that  were  observed  in  previous  studies  [4,26-28],  e.g.,  ethyl¬ 
benzene,  ethenyl-benzene  (styrene),  benzyl  alcohol,  benzoic  acid, 
dibenzyl,  and  methylbutenylbenzene.  Ethyl-benzene  was  identified 
but  not  quantified  as  it  coeluted  with  2-isopropyl-3,3-dimethyl- 
oxetane.  However,  species  such  as  benzyl  alcohol,  and  benzoic  acid, 
and  other  hydrocarbons  with  boiling  points  above  450  K,  were 
not  identified  due  to  chromatographic  limitations.  To  determine 
if  benzaldehyde,  benzene,  and  phenol  were  the  principal  species 
resulting  from  the  oxidation  of  the  toluene  and  the  extent  of  pos¬ 
sible  ring  rupture,  a  carbon  “ring  balance”  was  determined.  The 
“ring  balance”  ranged  from  110%  to  95%  at  the  start  of  the  NTC 
region.  The  apparent  ring  creation  was  an  artifact  of  the  previ¬ 
ously  mentioned  variation  in  the  toluene  concentration  between 
the  fuel  calibration  step  and  the  experiments.  The  completeness 
of  the  “ring  balance”  suggested  that  the  aromatic  ring  structure 


remained  intact  throughout  the  low  temperature  regime.  In  ex¬ 
periments  at  slightly  higher  temperatures,  Dagaut  et  al.  [29]  sug¬ 
gested  that  a  key  decomposition  pathway  for  phenyl  radicals  was 
through  cyclopentadienyl,  which  will  form  1,3-cyclopentadiene.  No 
1,3-cyclopentadiene  was  identified  although  the  chromatographic 
technique  could  measure  this  compound,  further  suggesting  that 
the  ring  structure  remained  intact. 

4.  Discussion 

Several  significant  differences  were  observed  between  the  neat 
hydrocarbon  experiments  and  the  surrogate  mixture  experiment. 
First,  the  marked  increase  in  the  2-methyl-5-ethyl-THF  concen¬ 
tration  observed  for  neat  n-heptane  was  not  observed  for  the 
surrogate,  rather  the  surrogate  mixture  produced  a  nearly  steady 
increase  in  concentration  to  720  K,  Fig.  16.  This  suggested  that 
the  buildup  of  heptyl  (C7H15)  and  heptylhydroperoxy  (C7H14OOH) 
radicals  was  slower  with  the  surrogate  mixture.  The  fact  that  the 
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Fig.  19.  Alkene  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  0  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  (B,  i  )  Ethene;  (A,—  —  ) 
1-propene;  (A, . ,)  2-methyl-l-propene;  (#,^—  .)  cis-  and  trans-2-heptene;  and  (□-.  i )  trans-3-heptene.  Symbols  are  experimental  results;  lines  are  model  predic¬ 

tions. 


Temperature  (K) 


Fig.  20.  Alkene  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  0  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  (B,  ■  )  4,4-Dimethyl-  1-pentene; 

(A,—  —  )  4,4-dimethyl-2-pentene;  (A, . .)  2,4-dimethyl- 1-pentene;  (•,—  .)  2,4-dimethyl-2-pentene;  (Q*»  1 1)  2, 4, 4- trimethyl- 1-pentene;  and  (O.whwhmwhmi)  2,4,4- 

trimethyl-2-pentene.  Symbols  are  experimental  results;  lines  are  model  predictions. 


model  predicted  this  behavior  further  suggests  that  the  toluene 
addition  was  not  responsible  for  this  change.  Furthermore,  the 
relative  production  of  two  2-methyl-5-ethyl-THF  isomers  in  the 
mixture  was  more  complicated  than  for  the  case  of  neat  n-heptane. 
At  the  lowest  temperatures,  specifically  640  and  650  K,  the  trans¬ 
isomer  was  only  slightly  favored.  Between  660  and  700  K,  the  rel¬ 
ative  formation  of  the  trans-isomer  increased  to  1.6:1,  the  same  ra¬ 
tio  observed  during  the  neat  studies.  At  temperatures  above  700  K, 
the  preference  for  the  trans-isomer  decreased  to  a  1.3:1  ratio.  This 
shift  in  the  relative  formation  of  the  isomers  coincided  with  the 
start  of  the  NTC  region,  700  K,  in  the  neat  studies,  yet,  for  the 
mixture  the  change  in  the  relative  formation  of  the  isomers  did  not 
shift  to  lower  temperatures  as  the  start  of  NTC  region  had  shifted. 
Secondly,  the  relative  importance  of  several  of  the  intermediates 
changed  significantly  between  the  neat  and  surrogate  mixtures. 
To  allow  comparison  between  the  surrogate  mixture  and  the  neat 


compounds,  select  species  were  identified  that  are  only  products 
of  a  single  component.  For  each  identified  species,  the  concentra¬ 
tion  of  that  species  was  multiplied  by  its  number  of  carbon  atoms 
and  then  normalized  by  the  initial  number  of  carbon  atoms  from 
the  fuel  component  that  was  responsible  for  the  identified  species 
formation,  either  neat  or  in  the  surrogate  mixture.  The  compari¬ 
son  of  this  value,  herewith  referred  to  as  “carbon  fraction,”  was 
made  between  the  neat  and  surrogate  mixtures  at  the  same  tem¬ 
peratures  during  the  experiments.  In  addition,  this  carbon  fraction 
was  also  calculated  and  compared  for  the  model  predictions  for 
the  neat  and  surrogate  mixtures  to  isolate  how  residence  time  and 
equivalence  ratio  differences  between  the  neat  and  surrogate  mix¬ 
ture  effected  the  observed  changes.  For  iso-octane  the  comparisons 
were  made  at  670  K  for  2,2,4,4-tetramethyl-THF,  2-isopropyl-3,3- 
dimethyl-oxetane,  2-methyl-propanal,  2,2-dimethyl-propanal,  2- 
methyl-l-propene,  2,2,4-trimethyl-l-pentene,  and  2,2,4-trimethyl- 
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Temperature  (K) 


Fig.  21.  Toluene  intermediate  species  profiles  as  a  function  of  reactor  temperature  for  GS4c-l  (5039  ppm,  </>  =  0.73,  65%  dilution,  225  ms,  and  8  atm):  (■)  Benzene;  (A ) 
benzaldehyde;  and  (♦ )  phenol. 


2-pentene.  For  n-heptane  the  comparisons  were  made  at  700  K  for 
2-methyl-5-ethyl-THF,  2-propyl-THF,  cis-  and  trans-2-heptene,  and 
trans-3-heptene. 

For  the  iso-octane  experiment,  the  measured  carbon  fraction 
of  2,2,4,4-tetramethyl-THF,  2-isopropyl-3,3-dimethyl-oxetane,  2- 
methyl-propanal,  and  2,2-dimethyl-propanal  increased  significantly 
by  29%,  31%,  70%,  and  45%,  respectively,  for  the  surrogate  mixture. 
The  experimental  carbon  fraction  for  the  alkenes  produced  from 
iso-octane  also  showed  a  sizable  increase,  2-methyl-l-propene, 
2,2,4-trimethyl-l-pentene,  and  2,2,4-trimethyl-2-pentene  increased 
48%,  16%,  and  16%,  respectively.  However,  when  the  increases  in 
the  experimental  carbon  fractions  were  compared  to  the  pre¬ 
dicted  carbon  fractions,  a  significant  difference  was  observed.  The 
model  predicted  more  modest  increases  in  carbon  fraction  of  16%, 
20%,  21%,  4%,  and  12%  for  2,2,4,4-tetramethyl-THF,  2-isopropyl-3,3- 
dimethyl-oxetane,  2-methyl-l-propene,  2,2,4-trimethyl-l-pentene, 
and  2,2,4-trimethyl-2-pentene,  respectively.  However,  the  forma¬ 
tion  of  the  iso-octane  derived  aldehydes  showed  the  most  signif¬ 
icant  differences  between  the  experimental  and  modeling  carbon 
fraction  as  the  model  only  predicted  a  9%  increase  in  2,2-dimethyl- 
propanal  and  no  significant  change  in  2-methyl-propanal. 

For  n-heptane,  the  experimental  carbon  fraction  for  2-methyl- 
5-ethyl-THF  and  2-propyl-THF  decreased  11%  and  22%,  and  cis-  and 
trans-2-heptene  and  trans-3-heptene  increased  42%  and  44%,  re¬ 
spectively.  Due  to  the  low  concentrations  of  pentanal,  the  effect 
on  the  aldehydes  of  n-heptane  could  not  be  inferred.  However,  the 
model  predicted  carbon  fractions  indicated  that  2-methyl-5-ethyl- 
THF,  2-propyl-THF,  2-heptene,  and  3-heptene  should  all  increase, 
specifically,  7%,  13%,  51%,  and  60%,  respectively.  As  shown,  the 
most  significant  differences  in  the  measured  and  predicted  carbon 
fractions  occurred  with  the  aldehydes  formed  from  iso-octane  and 
the  ethers  formed  from  n-heptane.  This  suggests  that  the  addition 
of  1-pentene  and  toluene  may  effect  branched  and  linear  alka¬ 
nes  differently.  It  also  suggests  that  some  key  interactions  are  not 
accounted  for  in  the  mechanism.  Although  speculative,  these  inter¬ 
actions  may  include  co-oxidation  reactions  between  iso-octane  and 
n-heptane  as  suggested  by  Andrae  et  al.  [30]  or  an  effect  associated 
with  toluene  addition. 

As  noted  previously,  one  of  the  most  significant  observations 
with  this  study  was  the  characterization  of  the  reaction  of  toluene 
with  the  other  surrogate  components  and  the  formation  of  aro¬ 


matic  intermediates.  Several  pathways  have  been  proposed  to  ex¬ 
plain  toluene  oxidation  in  the  low  and  intermediate  temperature 
regimes.  Specifically,  Ciajola  et  al.  [31]  suggested  several  pathways 
to  form  benzaldehyde  and  benzene  during  their  oxidation  of  n- 
heptane/toluene  mixtures.  Vanhove  et  al.  [4]  proposed  pathways 
to  produce  ethyl-benzene  and  methylbutenylbenzene  and  other 
co-oxidation  reactions.  However,  neither  study  observed  phenol. 
Although  speculative,  the  following  pathways  could  explain  the 
formation  of  phenol. 

0  +  OH*>0OH, 

0  +  ROOH  0OH  +  RO. 

In  the  surrogate  mixture,  there  are  significant  concentrations  of  OH 
produced  from  aldehydes  which  could  easily  react  with  the  long 
lived  phenyl  radicals  (0).  This  termination  step  would  have  a  sig¬ 
nificant  influence  on  the  abstraction  reactions  from  both  alkanes 
and  alkenes  and  on  OH  addition  to  the  double  bonds  of  alkene, 
as  they  are  sensitive  to  OH  concentrations.  The  second  proposed 
reaction  results  in  the  formation  of  an  alkoxy  radical.  The  alkoxy 
radicals  are  unstable  and  rapidly  decompose  to  form  smaller,  stable 
oxygenated  species,  principally  aldehydes  and  ketones,  and  smaller 
alkyl  radicals.  These  smaller  alkyl  radicals  then  tend  to  undergo  fi- 
scission  to  form  even  smaller  alkyl  radicals  and  alkenes.  As  a  result, 
an  increase  in  aldehyde  formation  and  alkenes  should  be  observed 
with  the  addition  of  toluene.  The  experimental  results  did  show 
an  increase  in  aldehyde  formation  as  compared  to  model  predic¬ 
tions,  lending  some  credence  to  the  proposed  pathways.  Vanhove 
and  coworkers  have  suggested  that  interaction  between  the  long 
lived  radicals  formed  during  toluene  reaction  and  other  larger  rad¬ 
icals.  However,  the  interactions  between  these  larger  radicals  were 
not  observed  during  this  study.  Nevertheless,  this  study  clearly 
shows  that  toluene  actively  reacts  in  the  low  and  intermediate 
temperature  regimes  with  the  radical  pool  generated  by  the  other 
compounds. 

5.  Summary  and  conclusions 

The  oxidation  of  a  proposed  gasoline  surrogate  was  investigated 
to  examine  the  interaction  of  the  constituents  when  examining  the 
surrogate  mixture.  The  gasoline  surrogate  was  a  mixture  of  4.6%  1- 
pentene,  49.6%  iso-octane,  14%  n-heptane,  and  31.8%  toluene.  Each 
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of  the  constituents  was  examined  individually,  nominally  from  600 
to  800  K  at  a  pressure  of  8  atm  under  dilute  conditions.  Both, 
1-pentene  and  n-heptane  showed  a  very  strong  NTC  behavior, 
where  the  start  of  NTC  occurred  at  710  K  and  703  K,  respectively. 
Iso-octane  also  exhibited  NTC  behavior  which  started  at  665  K. 
Toluene  was  not  reactive  at  the  conditions  examined.  The  surrogate 
exhibited  a  strong  negative  temperature  coefficient  (NTC)  behavior 
which  started  at  693  K.  Comparing  the  reactivity  behavior  of  the 
surrogate  components  with  a  surrogate  mixture  suggests  that  the 
temperature  for  start  of  the  NTC  region  for  blends  should  be  bound 
by  the  start  of  the  NTC  region  for  the  individual  components  of 
that  blend. 

Detailed  speciation  was  completed  throughout  the  NTC  tem¬ 
perature  range  to  identify  and  quantify  the  major  intermedi¬ 
ates  species.  For  1-pentene,  the  results  showed  that  the  fuel 
reacts  through  both  alkane  and  alkene  type  pathways  produc¬ 
ing  significant  concentrations  of  formaldehyde,  butanal,  acetalde¬ 
hyde,  ethene,  2-propenal,  and  1-propene.  The  preferred  alkene 
type  pathway  was  through  HO2  addition  to  the  double  bond,  the 
“Waddington”  mechanism,  which  was  responsible  for  the  forma¬ 
tion  of  significant  quantities  of  butanal  and  formaldehyde.  The  pre¬ 
ferred  alkane  type  pathways  were  through  abstraction  of  the  allylic 
hydrogen  and  subsequent  decomposition.  For  n-heptane,  over  80% 
of  the  fuel  decomposed  to  intermediates  at  the  start  of  the  NTC 
region  and  produced  significant  concentrations  of  formaldehyde, 
acetaldehyde,  2-methyl-5-ethyl-THF,  ethene,  propanal,  and  butanal. 
For  iso-octane,  only  20%  of  the  fuel  decomposed  to  intermediates 
at  the  start  of  the  NTC  region  and  produced  2,2,4,4-tetramethyl- 
THF,  2-isopropyl-3,3-dimethyl-oxetan,  formaldehyde,  2-methyl-l- 
propene,  2-propanone,  and  2-methyl-propanal.  For  the  surrogate, 
the  carbon  balance  was  very  good  (>90%)  considering  the  com¬ 
plexity  of  the  mixture.  Most  of  the  intermediates  identified  were 
the  same  as  those  present  in  the  neat  studies.  Comparisons  be¬ 
tween  the  concentrations  of  the  intermediate  species  in  the  sur¬ 
rogate  and  neat  component  experiments  showed  that  iso-octane 
produced  significantly  more  aldehydes,  alkenes,  and  ethers  in  the 
mixture,  while  the  n-heptane  intermediates  did  not  change  signifi¬ 
cantly.  This  implies  that  the  addition  of  1-pentene  and  toluene  may 
affect  linear  and  branched  alkanes  differently.  During  the  oxidation 
of  the  surrogate,  over  14%  of  the  toluene  reacted  to  intermedi¬ 
ates  while  no  reactivity  was  observed  with  neat  toluene.  These 
intermediates  appeared  to  retain  the  aromatic  ring  structure  as 
no  linear  alkenes  indicative  of  ring  rupture  were  identified.  Four 
major  species  were  attributable  to  toluene  decomposition,  specifi¬ 
cally,  benzaldehyde,  benzene,  phenol,  and  ethyl-benzene.  Two  new 
pathways  have  been  suggested  to  account  for  the  formation  of  the 
phenol,  which  has  not  been  observed  previously  in  toluene  studies. 
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